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Abstract Modem hydrodynamical simulations offer nowadays a powerful means to trace 
the evolution of the X-ray properties of the intra-cluster medium (ICM) during the cos- 
mological history of the hierarchical build up of galaxy clusters. In this paper we review 
the current status of these simulations and how their predictions fare in reproducing the 
most recent X-ray observations of clusters. After briefly discussing the shortcomings of the 
self-similar model, based on assuming that gravity only drives the evolution of the ICM, 
we discuss how the processes of gas cooling and non-gravitational heating are expected 
to bring model predictions into better agreement with observational data. We then present 
results from the hydrodynamical simulations, performed by different groups, and how they 
compare with observational data. As terms of comparison, we use X-ray scaling relations 
between mass, luminosity, temperature and pressure, as well as the profiles of temperature 
and entropy. The results of this comparison can be summarised as follows: (a) simulations, 
which include gas cooling, star formation and supernova feedback, are generally successful 
in reproducing the X-ray properties of the ICM outside the core regions; (b) simulations 
generally fail in reproducing the observed "cool core" structure, in that they have serious 
difficulties in regulating overcooling, thereby producing steep negative central temperature 
profiles. This discrepancy calls for the need of introducing other physical processes, such as 
energy feedback from active galactic nuclei, which should compensate the radiative losses 
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of the gas with high density, low entropy and short cooling time, which is observed to reside 
in the innermost regions of galaxy clusters. 

Keywords Cosmology: numerical simulations • galaxies: clusters • hydrodynamics • 
X-ray: galaxies 



1 Introduction 

Clusters of galaxies form from the collapse of exceptionally high density perturbations hav- 
ing typical size of ^ 10 Mpc in a comoving frame. As such, they mark the transition be- 
tween two distinct regimes in the study of the formation of cosmic structures. The evolution 
of structures involving larger scales is mainly driven by the action of gravitational insta- 
bility of the dark matter (DM) density perturbations and, as such, it retains the memory of 
the initial conditions. On the other hand, galaxy-sized structures, which form from initial 
fluctuations on scales of ~ 1 Mpc, evolve under the combined action of gravity and of com- 
plex gas-dynamical and astrophysical processes. On such scales, gas cooling, star formation 
and the subsequent release of energy and metal feedback from supernovae (SN) and active 
galactic nuclei (AGN) have a deep impact on the observational properties of the diffuse gas 
and of the galaxy population. 

In this sense, clusters of galaxies can be used as invaluable cosmological tools and astro- 
physical laboratories (see Rosati et al. 2002 and Voit 2005 for reviews). These two aspects 
are clearly interconnected with each other. From the one hand, the evolution of the pop- 
ulation of galaxy clusters and their overall baryonic content provide in principle powerful 
constraints on cosmological parameters. On the other hand, for such constraints to be robust, 
one has to understand in detail the physical properties of the intra-cluster medium (ICM) 
and its interaction with the galaxy population. 

The sim plest rnodel to predict the properties of the ICM and their evolution has been 
proposed bv lKaiseil llQSfj). This model is based on the assumption that the evolution of the 
thermodynamical properties of the ICM is determined only by g ravity, with gas hea ted to the 
virial temperature of the hosting DM halos by accretion shocks ( iBvkov et alj2008l - Chapter 
7, this volume). Since gravitational interaction does not introduce any preferred scale, this 
model has been called "self-similar'Q. As we shall discuss, this model provides precise pre- 
dictions on the shape and evolution of scaling relations between X-ray luminosity, entropy, 
total and gas m ass, which have been tested against numerical hydrodynamical simulations 
(e.g.. lEke et al.i ri998; Brvan & Norman 1998). These predictions have been recognised for 
several years to be at variance with a number of observ ations. In particular, the observed re- 
lation between X-ray luminosity and temperature (e.g. lMarkevitcdll998l ; lAmaud & Evrard 
1 19991 : lOsmond & Pon man 200^. is steeper and the measured level of gas entropy higher 
than expected (e.g., Ponm an et al ] |2003l ; (pratt & Amaudll2005l) . especially for poor clusters 
and groups. This led to the concept that more complex physical processes, related to the 
heating from astrophysical sources of energy feedback, and radiative cooling of the gas in 
the central cluster regions, play a key role in determirung the properties of the diffuse hot 
baryons. 

' Strictly speaking, self-similarity also requires that no characteristic scales are present in the underlying 
cosmological model. This means that the Universe must obey the Einstein-de-Sitter expansion law and that 
the shape of the power spectrum of density perturbations is a featureless power law. In any case, the violation 
of self-similarity introduced by the standard cosmological model is negligible with respect to that related to 
the non-gravitational effects acting on the gas. 
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Although semi-analytical approaches (e.g.. lTozzi & Normaiill200i] : IVoitll2005l . and ref- 
erences therein) offer invaluable guidelines to this study, it is only with hydrodynamical 
simulations that one can capture the full complexity of the problem, so as to study in detail 
the existing interplay between cosmological evolution and the astrophysical processes. 

In the last years, ever improving code efficiency and supercomputing capabilities have 
opened the possibility to perform simulations over fairly large dynamical ranges, thus al- 
lowing to resolve scales of a few kiloparsecs (kpc), which are relevant for the formation 
of single galaxies, while capturing the global cosmological environment on scales of tens 
or hundreds of Megaparsecs (Mpc), which are relevant for the evolution of galaxy clus- 
ters. Starting fro m first attempts, in which on ly simplifi ed heating schemes were studied 
("e.g. jNavarro et al . 1995: Bialek et al. 2001; Borgani et al. 2002J, a number of groups have 



studied the effect of introducing also cooling (e.g. Katz & Whit3l 19931 : iLewis etal ] l200Cll : 



iMuanwong et al.ll200l ':'Dave et al."2002':'Tornatore et al. 2003), of more realistic sources of 
energ y feedback (e.g., Borgani et al. 2004: Kav et al. 2007: N agai et ai] |2007a: Si iacki et alj 
l2007h . of thermal conducti on (e.g..lDolag et al.ll2004 ). and of non-t hermal pressure suppo rt 
from magnetic fields (e.g.. lDolag et al.ll2001 ) and cosmic rays (e.g. JPfrommer et"aDl2007h . 

In this paper, we will review the recent advancement performed in this field of computa- 
tional cosmology and critically discuss the comparison between simulation predictions and 
observations, by restricti ng the discussion to the thermal effects. As such, this paper com- 
plements the reviews bv lBorganieTai]|2008l - Chapter 18, this v olume, which reviews the 
study of the ICM chemical enrichment and by Dolag et al. 2008b - Chapter 15, this volume, 
which reviews the stu dy of the non-therm al properties of the ICM from simulations. We 
refer to the reviews bv lDolag et aDl2008al - Chapter 12, this vol ume for a description of the 
techniques of numerical simulations and bv lKaastra etai]|2008l - Chapter 9, this volume for 
an overview of the observed thermal properties of the ICM. 

The scheme of the presentation is as follows. In Sect. [2] we briefly discuss the self- 
similar model of the ICM and how the action of non-gravitational heating and cooling are 
expected to alter the predictions of this model. Sect. [3] and |4] overview the results obtained 
on the comparison between observed and simulated scaling relations and profiles of X-ray 
observable quantities, respectively. In Sect.[5]we summarise and critically discuss the results 
presented. 



2 Modelling the ICM 

2. 1 The self-similar scaling 

The simplest model to predict the observable properties of the ICM is based on the as- 
sumption that gravity only determines the thermodynamical properties of the hot diffuse gas 
Ikaiser 1986 ). Since gravity does not have preferred scales, we expect clusters of different 
sizes to be the scaled version of each other. This is the reason why this model has been called 
self- similar. 

If, at redshift z, we define M^^ to be the mass contained within the radius r^^ , encom- 
passing a mean density times the critical density Pc{z), then M^^ °c Pc(z)4cr^^. The 

critical density of the universe scales with redshift as Pc(z) = PcoE^{z), where E{z) is given 
by 

E{z)= [{i+z)'n,^ + {i+zfni^+nA]''\ (i) 
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where Q^n and are the density parameters associated to the non-relativistic matter and to 
the cosmological constant, respectively, Q]^ = I— — 12^ and we neglect any contribution 
from relativistic species. 

Therefore, the cluster size r^^ scales with z and M^^ as r^^ cx M^^E^'^I^ {z), so that, 
assuming hydrostatic equilibrium, cluster mass scales with temperature T as 

M^^T'I^E-\z). (2) 
If Pgas is the gas density, the corresponding X-ray luminosity is 

2 



-^-^ A(r)dV, (3) 
V \^lmp J 

where A{T) °c for pure thermal Bremsstrahlung emission. If gas accretes along with 
DM by gravitational instability during the formation of the cluster halo, then we expect that 
Pgas {r) °<: PDM{r), so that 

Lx^M^^p,T'/^^T^E{z). (4) 

Another useful qu antity characterising the thermodynamic al properties of the ICM is 
the entropy which, in X-ray studies of the ICM, is usually defined as 

where ka is the Boltzmann constant, /i the mean molecular weight (~ 0.58 for a plasma of 
primordial composition) and nip the proton mass. With the above definition, the quantity K 
is the constant of proportionality in the equation of state of an adiabatic mono-atomic gas, 

5/3 5/3 

P = Kpgis ■ Using the thermodynamic definition of specific entropy, s = cv^n{P/pgis ) (cy: 
heat capacity at constant volume), one obtains s = ke InK^^'^ + sq, where sq is a constant. 
Another quantity, often called "entropy" in the cluster literature, which we will also use in 
the following, is 

S^kBTn;^^\ (6) 

where «e is the electron number density. According to the self-similar model, this quantity, 
computed at a fixed overdensity A^, scales with temperature and redshift according to 

54c-r(l+z)-2. (7) 

As already mentioned in the introduction, a number of observational facts from X-ray 
data point a gainst th e simple self-similar picture. The s teeper slope of the Lx-T relation 
( lMarkevitchJ998 : Arnaud & Evrarj[l999lJosmond & Ponman 2004) . Lx °^ T" with a~3 
for clusters and possibly larger for groups, the excess entropy in poor clusters and groups 
JPonman et al]|2003l : IPratt & ArnaudI 12005 1; 'Piffaretti et al."2005') and the decreasing trend 
of the gas mass fraction in poorer systems (Lin et al. 2003; Sanderson et al. 2003) all point 
toward the presence of some mechanism which significantly affects the ICM thermodynam- 
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Grav. heat. SN feedb. S,i^50 kev cm^ 




Fig. 1 Profiles of reduced entropy, S/T, for non-radiative .simulations (from iBorgani et ai]l200lh . The left 
panel is for simulations including only gravitational heating, the central panel is for runs including SN feed- 
back as predicted by a semi-analytical model of galaxy formation and the right panel is with pre-heating wfith 
an entropy threshold at redshift z = 3. Solid, short-dashed and long-dashed curves are for a 3 keV cluster, 
for a 1 keV group and f or a 0.5 keV group, resp ectively. The dotted straight line in the left panel shows the 
analytical prediction bv lTozzi & Norman )200lh for the entropy profile associated to gravitational heating. 



2.2 Heating and cooling the ICM 

The first mechanism, tha t has been introduced to break the ICM self-similarity, is non- 
gravitational heating (e.g. lEvrard & Henrvlll99ll : lKaisedll99ll : lTozzi & Norman|[200lh . The 
idea is that by increasing the gas entropy with a given extra heating energy per gas particle 
iih prevents gas from sinking to the centre of DM halos, thereby reducing gas density and 
X-ray emissivity. This effect will be large for small systems, whose virial temperature is 
keT < E\^, while leaving rich clusters with Y^T ^> Zih almost unaffected. Therefore, we 
expect that the X-ray luminosity and gas content are relatively more suppressed in poorer 
systems, thus leading to a steepening of the Lx-T relation. 

The notion of non-gravitational heating has been first implemented in non-radiative (i.e. 
neglecting the effect of cooling) hydrodynamical simulati ons by either injecting entropy in 
an impulsive way at a given redshift ( Navarro et al. 1995i; lBialek et al.ll2001,) . or by adding 
energy in a redshift-modulated way , so as to mimic th e rate of SN explosions from an 
external model of galaxy formation l lBorgani et al.ll2002h . In Fig.[I] we show the different 
efficiency that different heating mechanisms have in breaking the self-similar behaviour of 
the entropy profiles in objects of different mass, ranging from a Virgo-like cluster to a poor 
galaxy group. According to the self-similar model, the profiles of reduced entropy, S/T , 
should be independent of the cluster mass. This is confirmed by the left panel, which also 
shows that these profiles have a slope consistent with that predicted by a model in which 
gas is shock heated by spherical accretion in a DM halo, under the effect of gravity only 
( iTozzi & Normanll200lh . The central panel shows instead the effect of adding energy from 
SN, whose rate is that predicted by a semi-analytical model of galaxy formation. In this 
case, which corresponds to a total heating energy of about 0.3 keV/particle, the effect of 
extra heating starts being visible, but only for the smaller system. It is only with the pre- 
heating scheme, based on imposing an entropy floor of 50 keV cm^, that self-similarity is 
clearly broken. While this heating scheme is effective in reproducing the observed Lx-T 
relation, it produc es large isentropic co res, a prediction which is at variance with respect to 
observations (e.g. jDonahue et al.ll2006h . 



6 




T (kev) log (i^gv) 

Fig. 2 Lef t panel: t he relation between entropy and temperature for gas having a fixed value of the cooling 
time (from lVoit &"B rvan 2001). The crosses with error bars are observational da ta on the entropy measured 
at one-tenth of the virial radius for clusters and groups bv iPonman et al J )1999 ). Right panel: a compari- 
son between observations (error bars with crosses; |PonmanetalJ<1999l) and simulations, including radiative 
cooling and star formation (numbers), for the entropy in the centr al regions of gala xy groups and clusters. 
The solid line shows the prediction of the self-similar model (from lDave et aU2002l) . 



Although it may look like a paradox, radiative cooling has been also suggested as a 
possible alternative to non-gravitational heating to increase the entropy lev el of the ICM 
and su ppressing the gas content in poor systems. As originally suggested bv lVoit & BrvanI 
coolin g provides a select ive removal of low-entropy gas from the hot X-ray emitting 
phase (see also Wu & XuelOOf). As a consequence, while the global entropy of the baryons 
decreases, the entropy of theX-i^y emitting gas increases. This is illustrated in the left panel 
of Fig.l2l(from lVoit & Brvanll200lh . In this plot, each of the two curves separates the upper 
portion of the entropy-temperature plane, where the gas has cooling time larger than the 
age of the system, from the lower portion, where gas with short cooling time resides. This 
implies that only gas having a relatively high entropy will be observed as X-ray emitting, 
while the low-entropy gas will be selectively removed by radiative cooling. The comparison 
with observational data of clusters and groups, also reported in this plot, suggests that their 
entropy level may well be the result of this removal of low-entropy gas operated by radiative 
cooling. This analytical prediction has been indeed confirmed by radiative hydr odynamical 
simula tions. The right panel of Fig. [2] shows the results of the simulations by lOave et al.l 
ll2002h on the temperature dependence of the central entropy of clusters and groups. Quite 
apparently, the entropy level in simulations is well above the prediction of the self-similar 
model, by a relative amount which increases with decreasing temperature, and in reasonable 
agreement with the observed entropy level of poor clusters and groups. 

Although cooling may look like an attractive solution, it suffers from the drawback that 
a too large fraction of gas is converted into stars in the absence of a source of heating energy 
which regulates the cooling runaway. Indeed, while observations indicat e that only about 10 
per cent of the baryon content of a cluster is in the stellar phase (e.g.. lBalogh et al. l l200ll : 
iLin et an i2003l) . radiative simulations, like those shown in Fig. |2l convert into stars up to 
~ 50 per cent of the gas. 
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Another paradoxical consequence of cooling is that it increases the temperature of the 
hot X-ray e mitting gas at the centre of clusters. This is shown in the left panel of Fig.[3](from 
IXomatoree t al. 2003), which compares the temperature profiles for the non-radiative run of 
a Virgo-like cluster with a variety of radiative runs, based on different ways of supplying 
non-gravitational heating. The effect of introducing cooling is clearly that of steepening 
the temperature profiles in the core regions, while leaving it unchanged at larger radii. The 
reason for this is that cooling causes a lack of central pressure support. As a consequence, 
gas starts flowing in sub-sonically from more external regions, thereby being heated by 
adiabatic compression. As we shall discuss in Sect.|4l this feature of cooling makes it quite 
difficult to reproduce the structure of the cool cores observed in galaxy clusters. 

Steepening of the central temperature profiles and overcooling are two aspects of the 
same problem. In principle, the solution to this problem should be provided by a suitable 
scheme of gas heating which regulates star formation, while maintaining pressurised gas 
in the hot phase. The right panel of Fig. [3] compares the temperature-density phase dia- 
grams for gas particles lying in the central region of an SPH-simulated cluster, when using 
two different feedback efficiencies. The two simulations include cooling, star formation and 
feedback in the form of galactic winds pow ered by SN explosions, following the scheme in- 
troduced bv lSpringel & HemquistI h003ih . The upper (green) and the lower (red) clouds of 
high temperature particles correspond to a wind velocity of 500 kms^' and of 1000 kms^' , 
respectively. This plot illustrates another paradoxical effect: in the same way that cooling 
causes an increase of the temperature of the hot phase, supplying energy with an efficient 
feedback causes a decrease of the temperature. The reason for this is that extra energy com- 
pensates radiative losses, thereby maintaining the pressure support for gas which would oth- 
erwise have a very short cooling time, thereby allowing it to survive on a lower adiabat. It is 
also worth remind i ng that cooling effi ciency increases with the numerical resolution (e.g., 
ISalogh et al.l200ll ; lBorgani et al.l2009) . Therefore, for a feedback mechanism to work prop- 
erly, it should be able to stabilise the cooling efficiency in a way which is independent of 
resolution. 

In the light of these results, it is clear that the observed lack of self-similarity in the 
X-ray properties of clusters cannot be simply explained on the grounds of a single effect. 
The emerging picture is that the action of cooling and of feedback energy, e.g. associated 
to SN explosions and AGN, should combine in a self-regulated way. As we shall discuss 
in the following sections, hydrodynamical simulations of galaxy clusters in a cosmological 
context demonstrate that achieving this heating/cooling balance is not easy and represents 
nowadays one of the most challenging tasks in the numerical study of clusters. 

As an example, we show in Fig. |4] how the gas density of a simulated cluster changes, 
both at z = 2 (left panels) and at z = (right panels), when coolin g and star formatio n 
are combined with different forms of non-gravitational heating (from lBorgani et al ] |200^ . 
The comparison of the top and central panels shows the effect of increasing the kinetic en- 
ergy carried by galactic outflows by a factor of six. The stronger winds are quite efficient 
in stopping star formation in the small halos, which are washed out, and make the larger 
ones slightly puffier, while preserving the general structure of the cosmic web surrounding 
the Lagrangian cluster region. Comparing the top and the bottom panels shows instead the 
effect of adding to galactic winds also the effect of an entropy floor. Although the energy 
budget of the feedback schemes of the central and bottom panels are quite comparable, the 
effect of the gas distribution is radically different. Imposing an entropy floor at z = 3 with 
an impulsive heating generates a much smoother gas density distribution, both at z = 2 and 
at z = 0. In this case, the filamentary structure of the gas distribution is completely erased, 
while only the largest halos are able to retain part of their gas content. This demonstrates 
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Fig. 3 Left panel: temperature profiles from hydrodynamical simulations of a ~ 3 keV galaxy cluster In all 
panels the dotted and the solid curves correspond to a non radiative run and t o a run including co oling and 
star formation. The other curves are for different recipes of gas heating (from lTomatore et al.|[200 3). Right 
panel: the relation between temperature and overdensity for gas particles within 0. lr2oo for SPH simulations 
of a cluster of mass ~ lO'^/i^'M ;,. Upper (green) points are for a run which includes feedback through 
galactic winds with a velocity of 500 kms"', while the lower (red) points are for a ran based on assum- 
ing stronger winds (sw), with a twice as large velocity. Both runs include a model of chemical enrichment 
jBorgani et al. 2008 - Chapter 18, this volume) which assumes an Initial Mass Function for star formation by 
lArimoto & Yoshik 19871 (AY). The points in the bottom right comer are star-forming gas particles. 



that a fixed amount of energy feedback can provide largely different results on the ICM ther- 
modynamical properties, depending on the epoch and on the density at which it is released. 



3 Scaling relations 

So far, we have qualitatively discussed how simple models of pre-heating and radiative 
cooling can reproduce the observed violation of self-similarity in the X-ray properties of 
galaxy clusters. In this and in the following sections we will focus the discussion on a more 
detailed comparison between simulation results and observational data, and on the implica- 
tions of this comparison on our current understanding of the feedback mechanisms which 
regulate star formation and the evolution of the galaxy population. As a starting point for the 
comparison between observed and simulated X-ray cluster properties, we describe how ob- 
servable quantities are computed from hydrodynamical simulations and how they compare 
to the analogous quantities derived from observational data. 

As for the X-ray luminosity, it is computed by summing the contributions to the emis- 
sivity, e,-, carried by all the gas elements (particles in a SPH run and cells in an Eulerian 
grid-based run), Lx = where the sum extends over all the gas elements within the 
region where Lx is computed. The contribution from the /-th gas element is usually written 
as 

£/ = «e./«H,;A {Ti,Zi)AVi , (8) 

where Wg,; and «h,i are the number densities of electrons and of hydrogen atoms, respectively, 
associated to the i-th gas element of given density p,-, temperature T, and metallicity Z,-. 
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Fig. 4 The maps of the gas density for simulations of a Viigo-hke cluster at ; = 2 and z. = (left and 
right panels, respectively), including cooHng, star formation and different forms of non-gravitational heating. 
Upper panels are for a run which includes galactic winds with a velocity of about 340 km s^', the central 
panels is for galactic winds with velocity of about 830kms"' and the bottom panels is for galactic winds 
as in the top panels, but also adding a pre-heating with an entropy threshold of 100 keV cm^ at z = 3. At 
z = the size of th e box is of 11.7/z"'Mpc, while at z = 2 is corresponds to 17.5/i^'Mpc comoving (from 
iBorgani et alj2005h . 
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Furthermore, AjT.Z) is th e temperature- and metallicity-dependent coohng function (e.g., 
[Sutherland & Dopital 19931 ) computed within a given energy band, while dV/ = ffi, /p, is the 
volume of the i-th gas element, having mass m, . 

As for the temperature, different proxies to its X-ray observational definition have been 
proposed in the literature, which differ from each other in the expression for the weight 
assigned to each gas element. In general, the ICM temperature can be written as 

T = ^ , (9) 

/ 

where 7] is the temperature of the i-gas element, which contributes with the weight w, . The 
mass-weighted definition of temperature, r^w, is recovered for w/ = m, (m,-: mass of the 
!-th gas element), which also coincides with the electron temperature for a fully ionised 
plasma. A more observation-oriented estimate of the ICM temperatur e is provided by th e 
emission-weighted definition, Tgw, which is obtained for w, = e, (e.g.. lEvrard et al.lll996h . 
The idea underlying this definition is that each gas element should contribute to the overall 
spe ctrum ac cording to its emissivity. 

iMazzotta et al, (.2004) pointed out that the thermal complexity of the ICM is such that 
the overall spectrum is given by the superposition of several single-temperature spectra, 
each one associated to one thermal phase. In principle, the superposition of several single- 
temperature spectra cannot be described by a single-temperature spectrum. However, when 
fitting it to a single-temperature model in a typical finite energy band, where X-ray tele- 
scopes are sensitive, the cooler gas phases are relatively more important in providing the 
high-energy cut-off of the spectrum and, therefore, in deter mining the temperature result- 
ing from the spectral fit. In order to account for this effect, iMazzotta et al. I ( |2004|) mtro- 
duced a spectroscopic-like temperature, T^\, which is recovered from Eq. [9]by using the 
weight Wi = PimiT"^^!'^ . By using a = 0.75, this expression for Tsi was shown to reproduce 
within few percent the temperature obtained from the spectroscopic fit, at least for clusters 
with temperature above 2-3 keV. More complex fitting expressions have been provided by 
IVikhlininl jloO^ , who generalised the spectroscopic-like temperature to the cases of lower 
temperature and arbitrary metallicity. 



3.1 The luminosity-temperature relation 

The Lx-T relation represented the first observational evidence against the self-similar model. 
This relation has been shown by several independent analyses to have a slope, Lx °<: T", with 
a ~ 3 for r > 2 keV (e.g.. White et al. 1997), with indications for a flattening to a < 2.5 for 
the most massive systems t Allen & Fabianll 19981) . The scatter in this relation is largely con- 
tributed by the coo l-core emission, so that it signific antly decreases when exc ising the cores 
l lMarkevitctJll998h or removing cool-core systems ( lAmaud & Evrardll 19991) . A change of 
behaviour is also observed at the scales of groups, T <2 keV, which generally displays a 

very large scatter (Osmond & Ponrnan 2004). 

Hydrodynamical simulations bv lBialeketaTl l l200lb and bv lBorgani et"ai] | |2002l) demon- 
strated that simple pre-heating models, based on the injection of entropy at relatively high 
redshift, can reproduce the observed slope of the Lx-T relation. JDave et al. (2002) showed 
that a similar result can also be achiev ed with simulatio n s incl uding cooling only, the price 
to be paid being a large overcooling. iMuanwong et al.l ( l2002h and iTomatore et alj ( l2003h 
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Fig. 5 Left panel: The relation between bolometric X-ray luminosity and emission-weighted temperature for 
simulations by Borgani et al. (2004), including cooling, star formation and feedback in the form of gal actic 
winds (green circles), compared with observational data for clusters I Amaud & Evrard 1999; Marl^e vitci] 
[l998) and groups (Helsdon & Ponman 2000). Right panel: the relation between X-ray luminosity, estimated 
outside the core regions, and the spectroscopic-like temperatu re, for simulations which include cooling, star 
formation and a form of "targeted" feedback (from lKav et aUi2Q07i) . Observational data for clusters are the 
same as in the left panel. 



demonstrated that combining cooling with pre-heating models can eventually decrease the 
total amount of stars to an acceptable level, while still providing a slope of the Lx-T relation 
close to the observed one. 

In Fig.|5]we show the comparison between observations and simulations in which the 
non-gravitational gas heating assumes the energy budget made available by the SN explo- 
sions, whose rate is computed from the simulated star formation rate. T he left panel shows 
the results from the SPH GADGET simulations by iBorgani et al.l | |2004 ). These simulations 
included the same model of kinetic feedback us ed in the simulations shown in the top panel 
of Fig. |4] This feedback model was shown by ISpringel & HemquistI ( l2003bl ) to be quite 
successful in producing a cosmic star formation history similar to the observed one. Quite 
apparently, these simulations provided a reasonable relation at the scale of clusters, T > 3 
keV, while failing to produce slope and scatter at the scale of groups. 

The right panel show the results bv lKavetal.l ( l2007h . which these authors plotted only 
for systems with T > 2 keV. Kay et al. also used SPH simulations based on the GADGET 
code, but with a different feedback scheme. In this scheme, energy made available by SN 
explosions is assigned in a "targeted" way to suitably chosen gas particles, which surround 
the star-forming regions, so that their entropy is raised in such a way to prevent them from 
cooling. Therefore, while the amount of energy is self-consistently computed from star 
formation, the way in which it is thermalised to the gas is suitably tuned. These simulations 

J iredic t a too high normalisation of the Lx-T relation. This result is interpreted bv lKav et alj 
2007h as due to the fact that their simulations produce too low temperatures for clusters of 
a given mass, as a consequence of the incorrect cool-core structure. 

Besides the slope of the local Lx-T relation, also its evolution carries information about 
the thermodynamical history of the ICM. Thanks to the increasing statistics of distant clus- 
ters observed in the last years with Chandra and XMM-Newton, a number of authors anal- 
ysed this evolution out to the highest redshifts, z— 1.3, where clusters have been detected so 
far. These analyses generally indicate that the amplitude of the Lx-T relation has a positive 
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Fig. 6 Left p anel: A comparison of the evolution of the Lx-T relation for t he simulated clusters by 
iBorgani et at 1 2004) and the Chandra data analysed by lEttori et al] ^20043) (from lEttori et alj |2004a). The 
"plus" symbols are the relation at z = 0, with the dotted and the dashed lines providing the best-fit relations 
when the slope is kept fixed at 2 and is left free, respectively. The squares, the triangles and the crosses 
show the simulation results at z = 0.5, 0.7 and 1, respectively. The small and the big dots show the observa- 
tional data at 0.5 < z < 0.8 and at z > 0.8, respectively. The lower panel shows the ratio between the mea- 
sured luminosities at z = ("plus" symbols) and at z = 1 (crosses) and the corresponding best-fitting power 
laws. Right panel : norm alisation of the Lx-T relation as a function of redshift, for the simulated clusters by 
iMuanwong et alj )2006l) in the radiative ntn with no feedback (solid line), in a run including an impulsive 
pre— heating at z = 4.5 (dashed fine), and in a run including a form of "targeted" feedback (dot-dashed fine). 



evolu tion out toz~ 0.5-0.6 (e.g. lAmaud et al.ll2002l : lLumb et al.ll2004l ; lKotov & Vikhliniii 



Ettori et al 



2005h, with hints for a po ssible inv ersion of this t rend at higher redshift (e.g 

2004b|; lMaugh an et al]|2006 : .Branchesi et alj|2007l) . In the attempt of interpreting these re 
suits, [Voit (2005^ showed that radiative cooling, combined with a modest amount of pre- 
heating with extra entropy, predicts an evolution of the Lx-T relation slower than that of the 
self-similar model, also with an inversion of the trend at high redshift. Although in quali- 
tative agreement with observations, these models have ho wever difficulties in fo llowing the 
observed positive evolution at z < 0.5 (e.g., see Fig. 14 bv lMaughan et al.ll2006h . 



As for the comparison with numerical predictions. lEttori et al. Ill2004ah analysed the evo 



lution of the Lx-T relation from the radiative simulation by iBorgani et al.l ( 12004') . which 
includes the effect of SN feedback in the form of galactic winds. As a result, they found that 
the normalisation of this relation for the s imulated clusters at high redshift is higher than 
for real clusters (see left panel of Fig. |6]l. iMuanwong et al.l ( 12006 ) analysed three sets of 
simulated clusters, based on radiative cooling only, on gas pre-heat ing at high redshift and 
on the same "targeted" SN feedback model used bv^ Kav et al.l ( 12007 '). As shown in the right 
panel of Fig. [6] they found that these three models predicts rather different evolutions of the 
Lx-T relation, thus confirming it to be a sensitive test for the thermodynamical history of 
the ICM. However, so far none of the numerical models proposed is able to account for the 
observational indication for an inversion of the Lx-T evolution at z ~ 0.5. This represents 
still an open issue, whose implications for a realistic modelling of the ICM physics are still 
to be understood. 
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3.2 The mass-temperature relation 

The relation between total collapsed mass and temperature has received much consider- 
ation both from the observational and the theoretical side, in view of its app lication for 
the use of gal axy clusters as tools to measure cosmological parameters (e.g., IVoitll200^ : 
lBorganill2003) . The relation between ICM temperature and total mass should be primar- 
ily dictated by the condition of hydrostatic equilibrium. For this reason, the expectation is 
that this relation should have a rather small scatter and be insensitive to the details of the 
heating/cooling processes. For a spherically symm etric system, t he condition of hydrostatic 
equilibrium translates into the mass estimator (e.g. lSarazinlll988h 



M{< r)= ^ ^ 



dlnpga,s(r) ^ dlnr(r) 



d In ?- d In r 



(10) 



Here M{< r) is the total mass within the cluster-centric distance r, while T(r) is the tem- 
perature measured at r. As shown in the previous section, the processes of heating/cooling 
are expected to modify the gas thermodynamics only in the central cluster regions, while the 
bulk of the ICM is dominated by gravitational processes. This implies that, in principle, the 
total mass estimate from Eq. [TO] should be rather stable. However, since the X-ray spectro- 
scopic temperature is sensitive to the thermal complexity of the ICM in the central regions 
(e.g., Mazzotta et al. 2004), a change of temperature and density profile in these regions may 
translate into a sizable effect in the mass-temperature relation provided by Eg. 1101 

Eq. llOlhas been often applied in the literature by modelling the gas density profile with 
a jS-model jCavaliere & Fusco-Femi ano 1976), 



-"gas ■ 



where rc is the core radius, and assuming either isothermal gas or a polytropic gas, p„ 
r^^', to account for the presence of temperature gradients. In this case, Eq.[TO]can be recast 
in the form 



M(<r) ^ l.lxlO'^jSr^^^^ ^T^- (12) 

^ ' ^' keV A-iMpc l-F(r/rc)2 



Ba sed on the analysis of non-radiative cluster simulations , I Schindlerl ( Il996h and lEvrard et all 
( Il996h argued that the X-ray temperature provides a rather precise determination of the clus - 
ter mass, with an intrinsic scatter of only about 15 per cent at rgoo.'Finoguenov et al.' ('2001') 
appli ed Eq.[T2]to ROSAT imaging and ASCA spectroscopic data (see also Nevalainen et al. 
l2000h . They found that the resulting M -T relation has a n ormalisation about 40 per cent 
lower than that from the simulations bvl Evrardet"an ( ll996l) . Although introducing the ef- 
fect of cooling, star formation and SN feedback provides a 20 per cent lower normalisa- 
tion, tiiis_wasjiorj^ enough^ with observations. Independent analy- 
ses dMuanwong et al.|[2002l ; Borgani et al.ll2004l) showed that applying to simulated clusters 



Eq. [12] which is used to estimate masses of real clusters, leads to a mass underestimate of 
about 20 per cent. This bias in the mass estimate is enough to bring the simulated and the 
observed M-T relations into reasonable agreement. These analyses were still based on the 
em ission-weighted te mperature in the simulation analysis. 

iRasia et al.l ( l2005l) showed that using the spectroscopic-like definition, T^\, leads to a 
mass underestimate of up to ~ 30 per cent with respect to the true cluster mass. This result is 
shown in Fig. [7] The left panel reports the relation between for simulat ed clusters and the 
true total cluster mass, also compared with the observed M-T relation bv lFinoguenov et aTl 
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Fig. 7 The relation between mass and spectroscopic-like temperature within rsoo- Red circles and green tri- 
angles are for simulations, which include cooling, star for mation and feed back from galactic winds, while 
squares with en'or bars are the observational data bv lFinoguenov et al.U200ll) . Left panel: M500 exactly com- 
puted by summing the mass of all the particles within rjoQ. Right panel: M500 estimated as in the observatio nal 
data, by using the equation of hydrostatic equilibrium for a polytropic /3-model. From lRasia et alji2005l) . 



( l200lh . Quite apparently, the relation from simulations lies well above the observational 
one. However, this difference is much reduced in the right panel, where the masses of the 
simulated clusters are computed by applying Eq.[T2] The reason for this difference between 
"true" and "recovered" masses is partly du e to the v iolation of hydrostatic equilibrium, as- 
sociated to subsonic gas bulk motions (e.g.. lRasia et al . 2004; Nagai et al. 2007b) and partly 
to the poor fit provided by the jS-model (e.g., Ascasibar et al.ll2003l) when extended to large 
radii. It is also interesting to note that the mass estimator of Eq.[T2]also under-predicts the 
intrinsic scatter of the M-T relation from the simulations. This is due to the fact that the 
equation of hydrostatic equilibrium imposes a strong correlation between ICM temperature 
and total mass. Any scatter is then associated to a cluster-by-cluster variation of the param- 
eters j6 and 7, which may not be fully representative of the diversity of the ICM structure 
among different objects. 

Using X MM-Newton and Chandra data, different authors (e.g., lAmaudetal1l2005l : 
Vikhlirun et al. 2005) applied the equation of hydrostatic equilibrium by avoiding the as- 
sumption of a simple beta-model for the gas density profile. As a result, the observed and 
the simulated M-T relations turned out to agree with each ot her, especially when simulated 
cluster s are analysed in the same way as real clusters (e.g.. lRasia et al ] l2006l : lNagaietalJ 
2007bh. The left pan el of Fig.[8]shows the comparison between the observational results by 



Amaud et al. I boost) a nd sim ulations. Here the discrepancy with respect to the non-radiative 



runs by lEvrard et al. is all eviated when including the effect of star formation and 

SN feedback (iBorgani et al.l[2004h . In a similar way, the right panel shows the compari- 
son between the ob served M-T relation by [vikhlinin et al. ( 2005) and the simulations by 
iNagai et ail ( l2007ah who computed cluster masses by using the same method applied to the 
Chandra data. This plot further demonstrates the good level of agreement between simulated 
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(keV) 

Fig. 8 Left panel: the observed mass-temperature relation at rjoo for nearby clusters observed with XMM- 
Newton, compared with M-T relations from simulations (from Amaud et al. 2005). The solid line is the best 
fit to observations, while the dotted line is the best fit only to a subsample of hot clusters. The short-dashed 
line is the relation from the non-radiative sim ulations bvjEv rard et al. 1 1996), while the long-dashed line 
is the relation from the radiative simulations bv lBorgani et alj l2004). Right panel: the M-T relation at rsoo 
from the radiative simulations by Nagai et al. 1 20 07al) : circles: z = 0; squares: z = 0.6, which include cooling, 
star-formation and an in efficient for m of SN feedback, compared to Chandra observations for a set of nearby 
relaxed clusters by Vik hlinin et alj (2005); stars with error bars. Masses of simulated clusters have been 
computed by using the same estimator, based on the hydrostatic equilibrium, applied to the Chandra data. 
Masses of both simulated and observed clusters have been rescaled with redshift according to the evolution 
predicted by the self-similar model. 

and observed M-T relation, once the 15-20 percent violation of hydrostatic equilibrium is 
taken into account. 



3.3 The mass-"pressure" relation 

To first approximation, the ICM can be represented by a smooth gas distribution in pressure 
equilibrium within the cluster potential well. Therefore, the expectation is that gas pres- 
sure should be the quantity that is more directly correlated to the total collapsed mass. For 
this reason, any pressure-related observational quantity should provide a robust minimum- 
scatter proxy to the cluster mass. One such observable is the Comptonisation p arameter, 
measured through the Sunyaev-Zeldovich effect (SZE; e.g., ICarlstrom et"ani2002l . for a re- 
view), which is proportional to the ICM pressure integrated along the line of sight. Ob- 
servations of the SZ effect are now reaching a high enough quality for extended sets of 
clusters, to allow correl ating the SZ signal with X-ray observable quantities. For instance, 
iBonamente et alj ilOO f) analysed a set of 38 massive clusters in the redshift range 0. 14 < 
z < 0.89, for which both SZE imaging from the OVRO/BIMA interferometric array and 
Chandra X-ray observations are available. As a result, they found that the slope and the 
evolution of the scalings of the Comptonisation parameter with gas mass, total mass and 

X-ray temperature are all in agreement with the prediction of the self-simil ar model. 

In an attempt to provide an X-ray observable related to the pressure, iKravtsov et alj 

introduced the quantity Yx = M^-^^^T, defined by the product of the total gas mass 
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Fig. 9 Left panel: the relation between Yx = Mg^Tx and Msoo- Points with error bars are for Chandra obser- 
vational data (from lNagai et al.l2007j) . Masses of simulated clusters have been computed by using the same 
estimator, based on the hydrostatic equilibrium, applied to the Chandra data. Masses of both simulated and 
observed clusters have been rescaled with redshift according to the evolution predicted by the self-similar 
model. Right panel: the relation between Yx and the X-ray luminos ity, both estima ted within r^oo^ lor a set of 
clusters at 0.1 < z < 1.3 extracted from the Chandra archive (from MaughanlllOOTl) . Luminosities have been 
rescaled according to the redshift dependence expected from self-similar evolution. 

times the temperature, both measured within a given aperture. With this definition, Yx repre- 
sents the X-ray counterpart of the Compton-y parameter, measured from the SZ effect. By 
computing this quantity for a set of simulated clusters, Kravtsov et al. (2006) showed that 
Yx has a very tight correlation with the cluster mass, with a remarkably small scatter of only 
8 per cent. The application to Chandra observations of a set of nearby relaxed clusters also 
shows that this relation has a comparably small scatter. Again, once masses of simulated and 
observed clusters are computed by applying the same hydrostatic estimator, th e normalisa- 
tion of the Yx-M relation from models and data closely agree with each other jNagai et"al] 
l2007al : see left panel of Fig. |9ll. Since mass is determined in this case from temperature, Yx 
and M are not independent quantities and, therefore, the scatter in their scaling relation may 
be underestimated. 

lMaughai3 l l2007h computed Yx for an extended set of clusters extracted from the Chandra 
archive, in the redshift range 0.1 < z < 1.3. The results of his analysis, shown in the right 
panel of Fig.|9l indicate that Yx is also tightly correlated with the X-ray luminosity through a 
relation which evolves in a self-similar way. Again, since the gas mass is obtained from the 
X-ray surface brightness, Yx and Lx are not independent quantities, thus possibly leading to 
an underestimate of the intrinsic scatter in their scaling relation. 



3.4 The gas mass fraction 

The measurement of the baryon mass fraction in nearby galaxy clusters has been recognised 
for several years to be a powerful method to measure the cosmological density parameter 
(e.g., lwhite et al.llI993l : iMotir et al.lil999l) , while i ts redshift evolution provides constr aints 
on the dark energy content of the Universe (e.g.. lAUen et alj|2002l ; lEttori et al.ll2003l and 
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Fig. 10 Left panel: the gas fraction (filled circles) and the baryon fraction (open circles) for cluster simula- 
tions which include cooling and star formation. Squares indicate runs where c ooling is switched off at z < 2, 
while triangles are for non-radiative simulations (froml Kravtsov et"ai]l2005h . Right panel: the gas fraction 
(diamonds) and the baryon fraction (circles) within rsoo for SPH simulations of clusters, which include dif- 
ferent gas physics. G: gravitational heating only; GV: G with a scheme for reduced SPH gas viscosity; FwW: 
radiative runs with feedback and weak galactic winds; FwWC: like FwW with also thermal conduction, with 
an efficiency of one-third of the Spitzer value; F: radiative runs with no winds; FsW: ra diative runs with 
strong winds. Blue, green and red symbols refer to redshifts 7, = 0. 0.7 and 1 (from Ettori et alj|20d3) . The 
shaded area indicates the range of values of the observed gas fraction tEttori & Fabiar^. 19991) . 



references therein). Since diffuse gas dominates the baryon budget of clusters, a precise 
measurement of the ICM total mass represents a fundamental step in the application of this 
cosmological test. While this method relies on the basic assumption that all clusters con- 
tain baryons in a cosmic proportion, a number of observ ational evidences show that the 
gas m ass fraction is smaller in lower temperature systems jLin et al ] |2003l :[ Sanderson et alJ 
I2OO3I) . This fact forces one to restrict the application to the most massive and relaxed sys- 
tems. In addition, since X-ray measurements of the gas mass fraction are generally available 
only out to a fraction of the cluster virial radius, the question then arises as to whether 
the gas fraction in these regions is representative of the cosmic value. Indeed, observa- 
tional evidence has been found for an increase of the gas mass fraction with radius (e.g., 
ICastillo-Morales & Schindierll2003h . 

In this respect, hydrodynamical simulations offer a way to check how the gas mass is 
distributed within individual clusters and as a function of the c luster mass, thus possibly 
providing a correctio n for such b iases. Indeed, ^llen et all ( |2004|) resorted to the set of SPH 
clusters simulated by lEke et ahl (1998.) to calibrate the correction factor that one needs to 
apply to extrapolate the baryon fraction computed at ^2500, which is the typical radius at 
which it is measured, to the cosmic value. However, since this set of simulations does not 
include the effects of cooling, star formation and feedback heating, it is not clear whether 
they provide a reliable description of the gas distribution within real clusters. 

iKravtsov et al.l ( |2005|) used high resolution simulations, based on an Eulerian code, for 
a set of clusters using both non-radiative and radiative physics. They found that including 
cooling and star formation has a substantial effect on the total baryon fraction in the central 
cluster regions, where it is even larger than the cosmic value. As shown in the left panel 
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of Fig. [To] at the virial radius the effect is actually inverted, with the baryon fraction of 
radiative runs lying below that of the non-radiative runs. In addition, they also compared 
results obtained from Eulerian and SPH codes and found small, but systematic, differences 
between the resulting baryon fractions. 

lEttori et al] feOOfit) performed a similar test, based on SPH simulations, but focusing on 
the effect of changing in a number of ways the description of the relevant physical processes, 
such as gas viscosity, feedback strength and thermal conduction. The results of their anal- 
ysis, which is shown in the right panel of Fig. (TO] demonstrated that the baryon fraction is 
generally stable but only at rather large radii, > rsoo- They also showed that changing the 
description of the relevant ICM physical processes changes the extrapolation of the baryon 
fraction from the central regions, relevant for X-ray measurements, while also slightly af- 
fecting the redshift evolution. 

On the one hand, these results show that simulations can be used as calibration instru- 
ments for cosmological applications of galaxy clusters. On the other hand, they also demon- 
strate that for this calibration to reach the precision required to constrain the dark energy 
content of the Universe, one needs to include in simulations the relevant physical processes 
which determine the ICM observational properties. 



4 Profiles of X-ray observables 

4. 1 The temperature profiles 

Already ASCA observations, despite their modest spatial resolution, have established that 
most of the clusters show significant departures from isothermality, with negative temper- 
ature gradient s characterised by a rem arkable degree of similarity, out to the largest radii 
sampled (e.g.. lMarkevitch et alj|l998|). Besides confirming the presence of these gradients. 



Beppo-SAX observations (e.g. lDe Grandi & M olendi 2002) showed that they do not extend 



down to the innermost cluster central regions, where instead an isothermal regime is ob- 
served, possibly followed by a decline of the temperature towards the centre, at least for re- 
laxed clusters. The much improved sensitivity of the Chandra satellite provides now a more 
detailed picture of the central temperature profiles (e.g., Vikhlinin et al. 2005; Baldi et al.l 
l2007h . At the same time, a number of analyses of XMM-Newton observatio ns now consis- 
tently show the presence of a negative gradient at radii > 0. lr2oo (e.g.. iPiffaretti et al I2OO5I : 
IPratt e t al. 2007, and references therein). Relaxed clusters are generally shown to have a 
smoothly declining profile toward the centre, reaching values which are about half of the 
overall virial cluster temperature in the innermost sampled regions, with non-relaxed clus- 
ters having, instead, a larger variety of temperature profiles. The emerging picture suggests 
that gas cooling is responsible for the decline of the temperature in the central regions, 
while some mechanism of energy feedback should be responsible for preventing overcool- 
ing, thereby suppressing the mass deposition rate and the resulting star formation. 

As for hydrodynamical simulations, they have shown to be generally rath er successful 
in reproducing the declining temperature profiles outside the core regions (e.g.. lLoken et al.l 
I2OO 2l; lRoncarelli et al .I2OO6I) , where gas cooling is relatively unimportant. On the other hand, 
as we have already discussed in Sect. |2.2| including gas cooling has the effect of steepening 
the T- profiles in the core regions, in clear disagreement with observations. The problem 
of the central temperature profile s in radiative sirnulations has been consistently found by 
several independent analyses (e.g. IValdaminill2003l ; iBorgani et al1l2004l ; iNagai et alj|2007af) 
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Fig. 11 Left panel: the temperature profiles in AMR simulations of clusters (thick curves) and in real clusters 
(thin curves). The simulation results are the average over 16 numerical clusters. The solid line is for runs with 
cooling and star formation, while the dashed line is for non-radiative simulations. The observational curves 
are for clusters with different temperatures (from Nagai et al. 20()7a). Right panel: Temperature profiles at 
z = 1 (top) and at z = (bottom) from a set of SPH simulated clusters including cooling, star-formation 
and a "targeted" scheme of SN feedback (from Kav e t alj|2007l) . Median and 10/90 percentiles are shown. 
Solid and dashed lines are for irregular and regular clusters, respectively. The vertical dashed line indicates 
the smallest scale which is numerically resolved. In the bottom panel, the heavy green line shows the profile 
from Chandra observations of cool core cl usters (V ikhlinin et aL2005i) . while the yellow shaded area is for a 
XMM-Newton sample of nearby clusters jPratt et all2d07l) . 



and is interpreted as due to the difficulty that currently implemented feedback schemes have 
in balancing the cooling runaway. 

As an example, we show in the left panel of Fig. [TT| the comparison betwee n simu- 
lated and observed temperature profiles, recently presented by 'Nagai et al.l ( l2007a h. which 
is based on a set of clusters simulated with an Adaptive Mesh Refinement (AMR) code. 
This plot clearly shows that the central profiles of simulated clusters are far steeper than the 
observed ones, by an amount which increases when cooling and star formation are turned 
on. Although this result is in qualitative agreement with other results based on SPH codes, 
an eye-ball comparison with the left panel of Fig. [3] shows a significant difference of the 
profiles in the central regions for the non-radiative runs. Indeed, while SPH simulations 
generally show a flattening at r < O.lryir, Eulerian simulations are instead characterised by 
continuously rising profiles. While this difference is reduced when cooling is turned on, it 
clearly calls for the need of performing detailed comparisons between d ifferent simulatio n 
codes, in the spirit of the Santa Barbara Cluster Comparison Project (e.g. lFrenk et aljl999h , 
and to understand in detail the reason for these differences. 

As discussed above, the steep temperature profiles predicted in the central regions wit- 
ness the presence of overcooling. Viceversa, the fact that a feedback mechanism is able 
to produce the cor rect temperature profiles does not guarantee in itself that overcooling is 
prevented. Indeed, lKavetalj|2007h found that their "targeted" scheme of SN feedback pro- 
duces temperature profiles which are in reasonable agreement with observations (see the 
right panel of Fig. [TTl i. However, even with this efficient feedback scheme, the resulting 
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Stellar fraction within clusters was still found to be too high, > 25 per cent, thus indicating 
the presence of a substantial overcooling in their simulations. 

Clearly, resolving the discrepancy between observed and simulated central temperature 
profiles requires that simulations are able to produce the correct thermal structure of the 
observed "cool cores". This means that a suitable feedback should compensate the radiative 
losses of the gas at the cluster centre, while keeping it at about ~ 1/3 of the virial tempera- 
ture. A number of analyses converge to indicate that AGN should represent the natural solu- 
tion to this problem. Considerable efforts have been spent to investigate how cooling can be 
self-consistently regulated by feedback from a central AGN in static cluster potentials (e.g. , 
IChurazov etaTltoOll : lOmma et al.l2004 iBrighenti & Mathewsll2006l : IStemberg et al.l2007l, 
and reference therein), while only quite recently these studies have beeri exten ded to clusters 
forming in a cosmological context (e.g. jHeinz et al.l2006l : ISijacki et al.l2007h . Although the 
results of these analyses are quite promising, we still lack for a detailed comparison between 
observational data and an extended set of cosmological simulations of clusters, convincingly 
showing that AGN feedback is able to provide the correct ICM thermal structure for objects 
spanning a wide range of masses, from poor groups to rich massive clusters. 



4.2 The entropy profiles 

As already mentioned in Sect. ID a convenient way of characterising the thermodynamical 
properties of the ICM is through the entropy, which, in X-ray cluster studies, is usually de- 
fined as iS = T jnj^ (see the review bv r^l2005l for a detailed discussion about the role 
of entropy in cluster studies). Since the current numerical description of the heating/cooling 
interplay in the central cluster regions is unable to reproduce the observed temperature struc- 
ture, there is no surprise that simulations have difficulties also in accounting for the observed 
entropy structure. However, while the ICM thermodynamics is sensitive to complex phys- 
ical processes in the core regions, one expects simulations to fare much better in the outer 
regions, say at r > 0.2r2oo, where the gas dynamics should be dominated by gravitational 
processes. This expectation is also supported by the agreement between the observed and 
the simulated slope of the temperature profiles outside cluster cores. 

If gravity were the only process at work, then the prediction of the self-similar model 
is that S °^ T. Therefore, by plotting profiles of the reduced entropy, 5/7", one expects them 
to fall on the top of each other for clusters of different temperature s. However, observa- 
tions revealed that this is not the case. Indications from ASCA data jPonman et aDl2003h 
showed that entropy profiles of poor clusters and groups have an amplitude which is higher 
than that expected for rich clusters from the above scaling arg ument. This result has been 
su bsequently confirmed by the XMM-Newton data analysed bv lPratt & Amaud ilOO^ and 
bv'Pif faretti et alj l l2005h . Quite remarkably, a relatively higher entropy for groups is found 
not only in the central regions, where it can arise as a consequence of the heating/cooling 
processes, but extends to all radii sampled by X-ray observations, out to < O.Srgoo. These 
analyses consistently indicate that S<xT", with a ~ 0.65, instead of a = 1 as expected from 
self-similar scaling. 

IPormian et al.1 1 20031) and lVoit et al.l ( l2003h interpreted this entropy excess in poor clus- 



ters and groups as the effect of entropy amplification generated by shocks from smoothed 
gas accretion. The underlying idea is the following. In the hierarchical scenario for structure 
formation, a galaxy group is expected to accrete from relatively smaller filaments and merg- 
ing sub-groups than a rich cluster does. Suppose now to heat the gas with a fixed amount of 
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Fig. 12 Left panel: the profile of reduced entropy, S/T, for a series of SPH GADGET runs of one cluster and 
three groups, including cooling, star formation and feedback from galactic winds. The top panel is for runs 
in which the winds have a velocity of abou t 360 kms~*, while the bottom panel is for an extreme wind 
model with a velocity of 830 kms"' (from Ifiorgani et ai]|2005t) . Th e straight dotted line marks the slope 
S <x R ^ '^^ which is the best-fitting to the observed entropy profiles (Pr att & Amaudll2005t iPiffaretti et al.l 
|2005|) . Right panel: the relation between entropy computed at rsoo for clusters and groups identified in AMR 
ENZO cosmological simulations (from I Younger & Brvanll2007i) . compared with the observational results 
bv iPon man et al. (2003); symbols with error bars. Shown are four different cases of entropy injection at 
z = 10: no pre-heating (open triangles), 78 keV cm~ (filled circles), 155 keV cm^ (stai's), 311 keV cm^ (filled 
triangles). The solid line is a power-law fit to the self-similar prediction from the simulations. 



specific energy (or entropy). Such a diffuse heating will be more effective to smooth the ac- 
cretion pattern of a group than that of a rich cluster, just as a consequence of the lower virial 
temperature of the structures falling into the former object. In this case, accretion shocks 
take place at a lower density and, therefore, are more efficient in g enerating entropy. While 
predictions from the semi-analytical approach bv lVoit et al.l ( |2003|) are in reasonable agree- 
ment with observational results, one may wonder whether these predictions are confirmed 
by full hydrodynamical simulations. 

To tackle this problem, iBorgani et all | |2005|) performed a series of SPH hydrodynami- 
cal simulations of four galaxy clusters and groups, with temperature in the range 0.5-3 keV, 
using both non-radiative and radiative runs, also exploring a variety of heating recipes. As 
a result, they found that galactic ejecta powered by an even extremely efficient SN feedback 
are not able to generate the observed level of entropy amplification. This result is shown in 
the left panel of Fig. [12] where the profiles of reduced entropy are plotted for the simulated 
structures in the case of moderate (upper panel) and very strong (bottom panel) galactic 
outflows. Borgani et al. also showed that only adding an entropy floor at z = 3 provides 
an efficient smoothing of the gas accretion pattern and, therefore, a substantial degree of 
entropy amplifications (see Fig.|4ll. Although this diffuse pre-heating provides an adequate 
entropy amplification at large radii, the price to pay is that the entropy level is substantially 
increased in the central regions. This is at variance with high-resolution Chandra measure- 
ments of low e ntropy gas in the inne rmost cluster regions, where it reaches values as low as 
~ 10 keV cm^ dOonahue et aklbOOd) . 

A similar analysis has also been performed by lYounger & BrvanI l l2007h . who used 
a set of AMR non-radiative cosmological simulations performed with the ENZO code 
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llO'Sheaetal.l l2004h. with gas entropy boosted at high redshift, z= 10. In keeping with 
previous analyses, they found that pre-heated simulations are generally able to reproduce 
the o bserved luminosity-te mperature and mass-temperature relations. However, differently 
from iBorgani et all ( l2005t) , even the most extreme pre-heating scheme does not provide 
an appreciable degree of entropy amplification. This result is shown in the right panel of 
Fig- El The deviations from self-similarity for the entropy compu ted at r^on are a lways too 
small to reconcile the simulations with the observational data by ^ Ponman et al.1 ( l2003h . A 
possible reason for the different results with respect to the analysis by Borgani et al. is due 
to the higher redshift of pre-heating, z = 10 instead of 3, used by Younger & Bryan. Fur- 
thermore, grid-based codes are known to produce entropy profiles that, in the central part of 
the halos, r < Q.2r2m, are flatter than those produced by SPH codes (e.g., Voit et al. 20(3). 
Therefore, the question arises as to whether different numerical schemes of hydrodynamics 
react in different ways to pre-heating. It would be highly recommendable that future simu- 
lation comparison projects will include tests of how different codes behave in the presence 
of simple schemes of non-gravitational heating. 

5 Summary 

In this paper we reviewed the current status of the comparison between cosmological hy- 
drodynamical simulations of the X-ray properties of the intra-cluster medium (ICM) and 
observations. We first presented the basic predictions of the self-similar model, based on the 
assumption that only gravitational processes drive the evolution of the ICM. We then showed 
how a number of observational facts are at variance with these predictions, with poor clusters 
and groups characterised by a relatively lower density and higher entropy of the gas. This 
calls for the need of introducing some extra physical processes, such as non-gravitational 
heating and radiative cooling, which are able to break the self-similarity between objects 
of different size. The results based on simulations, which include these effects, can be sum- 
marised as follows. 

(1) The observed scaling relation between X-ray luminosity and temperature can be repro- 
duced by simulations including cooling only, but at the expense of producing a too large 
fraction of cooled gas. Introducing ad-hoc schemes of entropy injection at high redshift 
(the so-called pre-heating) can eventually produce the correct Lx-T relation. However, no 
simulations have been presented so far, in which a good agreement with observations is 
achieved by using a feedback scheme in which the energy release and thermalisation from 
SN or AGN is self-consistently computed by the simulated star formation or accretion onto 
a cosmologically evolving population of black holes. 

(2) Simulations which include cooling, star formation and SN feedback, correctly reproduce 
the observed mass-temperature and "pressure'-temperature relations. Quite interestingly, 
this agreement is achieved once the masses of simulated clusters are estimated by apply- 
ing the same mass estimators, based on the assumption of hydrostatic equilibrium, which 
are applied to the analysis of real clusters. The violation of hydrostatic equilibrium, related 
to the non-thermal pressure support from subsonic gas motions, would otherwise lead to a 
~ 20 per cent overestimate of the normalisation of the above scaling relations for simulated 
clusters. 

(3) Both non-radiative and radiative runs naturally predict the observed negative gradients 
of the temperature profiles outside the cluster core regions, r > 0.2r2oo. This suggests that 
the ICM thermal structure in these regions is indeed dominated by the action of gravitational 
processes, such as heating from shocks associated to supersonic gas accretion. 
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(4) Introducing cooling has the effect of steepening the temperature profiles in the inner- 
most regions, as a consequence of the adiabatic compression of inflowing gas, caused by the 
lack of pressure support. Therefore, gas in the cores of simulated clusters generally lies in a 
high-temperature phase, which is very well separated from the cold (~ 10"* K) phase. This 
is at variance with the observed "cool core" structure of real clusters, which show instead 
the presence of a fair amount of gas, down to a limiting temperature of ~ 1 /3 of the virial 
temperature, which formally has a rather short cooling time. The presence of this gas makes 
the observed temperature profiles of relaxed clusters to peak at r 0.2r200i thereby decreas- 
ing by about a factor two in the innermost sampled regions. The discrepancy between the 
simulated and the observed "cool core" calls for the need of introducing in cluster simula- 
tions a suitable feedback mechanism which is able to compensate the radiative losses, keeps 
pressurised gas in the central regions, and suppresses the mass deposition rate. 

(5) Simulations are generally rather successful in reproducing the observed slope of the en- 
tropy profiles outside the core regions, S °^ r" with a '^\. This slope is also close to that 
predicted by models based on gravitational heating from spherical accretion, thus lending 
further support to the picture that gravity drives the evolution of the ICM outside the core 
regions. Quite intriguingly, however, the normalisation of the profiles out to the largest sam- 
pled radii, < r^Qo, is observed to have a milder temperature dependence than expected from 
the self-similar model, T" with a ~ 0.65 instead of a ~ 1. This entropy excess in poor 
systems can be reproduced in simulations only by adding a diffuse pre-heating mechanism, 
which smoothes the pattern of gas accretion, thereby leading to an amplification of entropy 
generation from shocks in poorer systems. However, this mechanism is able to provide the 
correct explanation only for clusters with temperature T <3 keV. Generating entropy ampli- 
fication in hotter systems would require an exceedingly large amount of pre-heating, which 
would generate too shallow entropy profiles. 

hi general, the above results demonstrate the capability of cosmological hydrodynam- 
ical simulations to predict the correct thermodynamical properties of the ICM outside the 
central regions of galaxy clusters. On the other hand, even simulations based on an effi- 
cient SN feedback fail in preventing overcooling and providing the correct description of 
the thermodynamical ICM properties in the central regions. The "cool core" failure of sim- 
ulations based on stellar feedback is also witnessed by the exceedingly massive and blue 
central galaxies that are generally produced (see also lBorgaiu et al.ll2008l - Chapter 18, this 
volume). 

The fact that cool cores are observed for systems spanning at least two orders of mag- 
nitude in mass, from ~ lO^^M^^j to ~ lO'^Mj,, suggests that only one self-regulated mech- 
anism should be mainly responsible for this, rather than the combination of several mech- 
anisms possibly acting over different time-scales. Although AGN are generally thought 
to naturally provide such a mechanism, a detailed comparison between observational data 
and an extended set of cluster simulations, including this kind of feedback is still lacking. 
Furthermore, understanding in detail the role of AGN in determining the ICM properties 
requires addressing in detail two major issues. Firstly, the accretion process onto the central 
black hole involves scales of the order of the parsec, while the X-ray observable effects 
involve scales of ~ 10-100 kpc. This requires understanding the cross-talk between two 
ranges of scales, which differ by at least four orders of magnitude. Secondly, the total energy 
budget available from AGN is orders of magnitude larger than that required to regulate cool- 
ing flows. Therefore, one needs ultimately to understand the channels for the thermalisation 
of the released energy and how this naturally takes place without resorting to any ad-hoc 
tuning. These problems definitely need to be addressed by simulations of the next genera- 



24 



tion, which will be aimed at understanding in detail the evolution of the cosmic baryons and 
the observational X-ray properties of galaxy clusters. 
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